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Abstract 

We present the results on non-perturbative quantum gravity effects 
related to extra dimensions which can be comparable, in some cases, 
with the SM contributions, e.g. in lepton-lepton or lepton-nucleon 
scattering. The case of cosmic neutrino gravitational interaction with 
atmospheric nucleons is considered in detail. 
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1 Gravi-Reggeon effects in multidimensional 
scattering amplidutes 

During last years there is a growing practical interest in models with compact 
extra spatial dimensions. Their compactification radius, R c , varies from 1 fm 
to 1 mm, depending on a number of extra dimensions d = D — 4 PQ. The 
models predict massive Kaluza-Klein (KK) excitations of the graviton and 
KK modes of the SM fields (provided the latter are allowed to propagate 
in higher dimensions). If D-dimensional space-time has a flat metric, the 
coupling of the massive graviton modes with the SM particle is very weak and 
is defined by the Newton constant Gn = where M Pt is the reduced 

Planck mass. Nevertheless, in the case when SM particles are confined to a 4- 
dimensional flat "brane" , summing up the KK graviton excitations results in 
a D-dimensional gravitational coupling Go ~ l/M D +d , with a fundamental 
Planck scale M D of order 1 TeV PQ. 

Let us first consider the SM in D -dimensional flat space-time, D > 4, 
without gravity. Due to extra spatial dimensions, an effective "transverse 
interaction region" becomes larger than in four dimensions. One manifesta- 
tion of this is a modification of the Froissart-Martin upper bound in a flat 
space-time with arbitrary D dimensions [2j: 



\/s being a collision energy. The "transverse radius" in JIJ is given by 
Ro(s) = N(D) Ins / y/to, where to denotes the nearest singularity in the t- 
channel, assumed non-zero, while N(D) is some integer depending on D. 
It is interesting to see, what happens with scattering amplitudes when we 
replace infinite extra dimensions by compact ones? 

In Ref. [2] the Froissart-Martin bound was generalized for scattering in 
D -dimensional space-time with compact extra dimensions. For one extra 
dimension with the compactification radius R c , the upper bound is of the 
form: 



where lmT(s,t) is the scattering amplitude, t is a momentum transfer (in 
D dimensions) and <$>(Rq/ R c , D) is a known function. At R c <C Ro(s) the 
equality (J2J results in [3] 
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while in the opposite limit, R c 3> -Ro( s ), the inequality (J2J) reproduces the 
upper bound (JTJ). 

Now let us allow for the gravity to come into play. As was argued in 
a number of papers in the Minkowski space-time with D dimensions 

(D > 4) the gravity becomes strong in a transplanckian region (s 3> Md), 
since an effective gravitational coupling, Gps, rises with energy 

In Refs |3] the eikonal representation for the scattering amplitude of the 
gravitons in the string theory was obtained: 

^ t ) = -«./^[^)-l], (4, 

where — I m x(s,&) is large at b < b\ = 2y/a^lns (a' g is a string 

tension). Thus, one get asymptotically 

ag(s) ^const(D) bf~ 2 (s). (5) 

Due to the absence of infrared divergences in the flat space-time with more 
than four dimensions, the gravitational cross section (J5J) appears to be finite 
and similar to the upper bound (JTJ). 

In what follows, we will first consider the scattering of two particles in 
the model with one compact extra dimensions (D = 5) in the transplanckian 
kinematical region: 

y/s > Md, s > -t, (6) 

t = —q]_ being four- dimensional momentum transfer. The generalization to 
D > 5 is straightforward and it will be done below. Thus, we start from the 
consideration of the scattering of bulk particles in four spatial dimensions, 
one of which is compactified with the large radius R c . 

In the eikonal approximation an elastic scattering amplitude in the trans- 
planckian kinematical region is given by the sum of reggeized gravitons 
in t-channel. So, we assume that both massless graviton and its KK massive 
excitations lie on linear Regge trajectories: 

a(t D ) = a(0) + a' g t D , (7) 

where tp denotes D-dimensional momentum transfer. Since the extra di- 
mension is compact with the radius R c , we come to splitting of the Regge 
trajectory (J7J) into a leading vacuum trajectory [H] 

a (t) = a grav (t) = 2 + a' g t (8) 
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and infinite sequence of secondary, "KK-charged" , gravireggeons: 



a'. 



a n (t) = 2--^n 2 + a' g t, n^l. (9) 



c 



The string theory implies that the slope of the gravireggeon trajectory is 
universal for all s, and a' g = 1/M S 2 , where M s is a string scale. 

If we assume that multidimensional theory at short distances is a string 
theory, than the scale Mp can be of order the fundamental string scale M s = 
(a') -1 / 2 . For instance, in the type I theory of open and closed strings one 
has 

, 2 \ 2/(2+d) 

M s = j M D , (10) 

where g s is a gauge coupling at the string scale. This relation leads to D- 
dimensional Planck mass a bit higher than the string scale (for g 2 /47r ~ 0.1). 

Thus, instead of taking a "bare" graviton exchange, we calculate a contri- 
bution from the Pomeron trajectory to which this KK graviton mode belongs: 

i + mp( -y)) am f±y"\ 

smira n (t) y \s ) 

The Born amplitude is, therefore, of the form 

r 7r i / s \ an{t) 

A B (s, t, n) = G N (2rrR c ) - cot -a n {t)\ a' g (3 2 n {t) I- . (12) 

In order to get an idea of possible t-dependence of Regge residues /? 2 (t), 
we consider scattering of D-dimensional gravitons. The corresponding am- 
plitude has been calculated in Refs. 0j: 



AB ^^^w^) {a '- s)<tD - (13) 



The expression fllH)l is valid in the region a' g \t\ < 1 in which it can be recast 
in the form: 

Af trm9 (s,t D ) ~^e^» (a' gS y'^, (14) 



where 7 ~ 0.58 is the Euler constant. 
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Thus, we have A(s,t) ~ exp(a' g ct), where c is of order of unity. Let us 
assume that Regge residues in (jT2j) have an analogous t-dependence: 

(? n (t) = (3 2 (0)e< bo{t - n2/R '\ (15) 

Since the coupling of all KK states to the SM fields is universal in ADD model, 
we expect that /3^(t) depends on n via to = t—n 2 /R 2 . Accounting for the fact 
that the product a' g bo appears only in a combination with a' g ln(s/so), we can 
neglect it in forthcoming calculations at large s and put /3 2 (t) ~ /3 2 (0). At 
t — > 0, n = expression (jTTj) should reproduce singular term Gjys/\t\ related 
with the massless graviton, that results in the relation 2/3 2 (0)/7tSq = 1. 

The expression for 5-dimensional eikonal amplitude looks like (k being 
the exchanged KK quantum number) 

7T 

A{s, t, k) = 2iR c s J d 2 b e iqi - h+ik * Jd<p[l- e ix{sM) ] , (16) 
with the eikonal given by 

oo 

i r , oo 

X (s, b,4>) = —l q±dq ± J (q ± b) — ^ e**A B ( 8 , -q 2 ± , n). (17) 

q c n=— oo 

The variable <fi runs the region — 7r ^ <p ^ ir. These inequalities imply that 
— oo ^ y ^ oo in the limit R c —>■ oo (flat extra dimension), y = R c <p being 
the 5-th component of the impact parameter. 

One can easily obtain from (|16|) that at k = and s < A/R 2 only modes 
with n = contribute and effectively x( s i b, (p) = Xo( s ^b), corresponding to 
n = contribution in the sum in Eq. (|17|). So, at low energy the scattering 
amplitude does not feel extra dimensions (the factor R c is trivial and is absent 
at proper normalization). 

Let us consider first the imaginary part of the eikonal. From equations 
(HZD, (H2D we obtain: 

a' r l 

Im X (s,M) = G^s^-yexp [-6 2 /4^(s)J 9 3 (v,q), (18) 

where 



R g (s) = Ja'(\n(s/s )+b ) (19) 
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is a gravitational slope. The quantity 9 3 in (jTHjl is one of Jacobi ^-functions ||: 

oo 

9 3 (v) = 9 3 (v, q) = 1 + 2 ^ cos(2™) q" 2 . (20) 



n=l 

In our case, it depends on variables 



V = 2? 
q = exp 



-R'M/K- (2i) 



The function 9${v,q) is well-defined for all (complex) v and all values of q 
such as \q\ < 1. It has a singularity at g — > 1 (see below). The ^-functions 
are often defined in terms of variable r: 

9(v)=9(v\r), (22) 

where 

q = e i7TT . (23) 

Let us define the ratio: 

(24) 



2i? 9 ( S ) 

(that is, q = exp(— l/4a 2 )). From the equality pQ 

/1TpV\ 1+2/d 

R c = 2-10 31 / d - 17 cm (25) 

we see that the compactification scale R~ x varies from 10 -3 eV for d = 2 to 
10 MeV for d = 6. Since i?^ 1 <C (2i2 s (s)) _1 even at ultra-high energies, we 
have a ^> 1 and, consequently, (1 - g) < 1. 

The behavior of Q${v,q) at q — > 1 can be derived by using unimodular 
transformation of ^3-function (known also as Jacobi imaginary transforma- 
tion) 0: 



- J = (-?r) 1/2 e i7r " 2/T ^ 3 (^|r). (26) 



Here (— ir) 1 ^ 2 has a principal value which lies in the right half-plane. In 
variable q, equality ()26|) looks like 



/ \ 1/2 00 



71= — CX) 
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The series in the RHS of (|27|) converges very quickly at q — > 1, contrary to 
original series (J2(J|) : 



9 3 {v, q) = 2a^\e-* 2a2 + ^ L-V™-^ + e -(a™ + *)V 



(28) 



n=l 



Notice that a 2 = -1/4 In g. 
From all said above, we get 



lmx(s, b, (f>) ~ G N s 



8m(s) 



cxp 



(29) 



The expression (I18J) is directly generalized for d extra dimensions (c/ ^ 1): 
lmx(s,&,0i,...0d 
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x exp 



8i^( s ) 
b 2 /4R 2 g (s)} Y[9 3 ( Vi ,q), 



(30) 



i=l 



where Vi = 4>i/2n. Correspondingly, we obtain 

a'RS d ' 2 



hn X {s, b, fa,... <j>i) ~ G N s % 9 R2 l d ^ 



x exp 



(b' + Ri^ + .-. + R^/ARiis) 



(31) 



We see from (13 1|) that the imaginary part of the eikonal decreases exponen- 
tially in variables b, (pi outside the region: 

b 2 + {R^f + ... {R c ct> d ) 2 < R 2 (s), (32) 

where 

R 2 (s)^AR 2 g (s)\n(s/M 2 D ) (33) 

at high s. 

Let tc = (t, —n\/R 2 c1 . . . , —n 2 d /R 2 ) be a bulk momentum transfer. Then 
we get the following expression for multidimensional scattering amplitude: 



A D (s,t,n 1 , ...,n d ) 



-2is R d c I d 2 b e iqi - b i u<pf- 



n 

8=1 



e ini4>i ^ e ix{s,b,4>l,---,(pd) _ ]_J 



(34) 



Correspondingly, the inelastic cross section in the space-time with d compact 
dimension is given by 



TT TT 

ag(s) = (2nR c ) d J d 2 b J dfa ■ ■ ■ J d(j> d 1 - e ^x(-M,-M _ (35) 

— TT —TT 

As was already shown, the imaginary part of the eikonal is negligibly small 
outside region (J3~2|) . That results in the estimates: 

r R 2 +d (s), i? c »i? (s) 

a? n {s) ~ const (D) x I (36) 

[ R 2 (s)R d c , R c <^R (s) 

which remind Eqs. (JIJ and Q obtained previously. As was mentioned 
above, Rq(s) <C R c for any reasonable s. So, the size of the compact ex- 
tra dimensions is irrelevant to the behavior of the inelastic cross section and 
a in( s ) ~ ( a ' g ) D ^ 21 (^ n S ) D 2 - Only at s — > 00, when the transverse interaction 
region Rq(s) becomes much larger than R c , we get cr^(s) ~ a' g R®~ A {\n.s) 2 . 

2 Scattering of the SM fields in the presence 
of compact extra dimensions 

Now we consider the case when the colliding particles are confined on the 
4- dimensional brane, while the exchange quanta (KK gravitons) are allowed 
to propagate in the bulk. Thus, the collisions of the SM particles take place 
in a two-dimensional impact parameter space. In Refs. HD1 the scattering 
of two SM particles was calculated in the eikonal approximation by summing 
up only "bare" KK gravitons. The massive graviton modes originated from 
the extra dimensions change four- dimensional propagator by 

nf+...n^0 - t + £ -L 

i=i R c 

Since a contribution from only non-reggeized KK excitations of the graviton 
has been taken into account in jUlEIl, the eikonal has no imaginary parts in 
such an approach. As was shown in ^U] , the D-dimensional brane amplitude 
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has a renormalized Born pole at t = and an infinite phase. Notice, series 
(|37|) diverges and needs renormalization at d ^ 2. 

In Ref. |TL it was shown that an amplitude of M — * N transition ob- 
served in four dimensions, Amn, is related to a corresponding D-dimensional 
amplitude Af. IN by the relation 

A MN = (2nR c ) d ^ M+N ^ A D MN (38) 

(in our case, M = N = 2). Amplitudes A^n have non-zero limit at R c — > 0, 
reproducing the usual 4-dimensional pseudoeuclidean case. Since the collid- 
ing particles are confined on the brane, their momenta lie in four-dimensional 
space. Therefore, the impact parameter belongs to the two-dimensional space 
and we have to put fa — 0, i — 1, . . . , d, in (j3*TJl . With taking account of this, 
the expression for four-dimensional eikonal amplitude (in the presence of d 
compact extra dimensions) looks like 

A(s, t) = 2is J d 2 b e iq±b [1 - e ix{s ' b) ] , (39) 

where x( s ; b) = x( s > b, 4>x — 0, . . . , <fid — 0). Taking into account (|3*T|) . we get 
the expression 



lmx(s,b) 



In 



xd/2-i M 2 \2M Dt 
x exp[-b 2 /AR 2 g ( S )], (40) 



where the relation Mp t = (2irR c ) d M^ d is used p. The detailed analysis of 
the real part of the eikonal will be given elsewhere. Here we only note that, 
contrary to ([411)1. the real part of the eikonal (with the massless graviton term 
subtracted) decreases as a power of the impact parameter at large b. 

The important features of the gravitational contribution to cross sections 
are its independence of types of colliding particles and a strong dependence 
on the collision energy . So, one can expect that at superplanckian energies 
gravity exchanges will dominate the SM electroweak interactions. That is 
why we now focus on leptonic and semileptonic collisions. 

Let us first consider e + e _ annihilation. Unfortunately, future linear col- 
liders will provide us with the c.m.s. energies \/~s around Md (0.5 -v- 2 TeV). 
In order to estimate af^ e numerically, we need to fix Regge free parameter s 
in (|40|) . Since Sq is related rather with a mass scale of exchange quanta than 
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with mass scales of colliding particles, we can treat the scattering amplitude 
of two graviton (J13j) instead of SM particle collision, and deduce that 1 

s = (of,)- 1 . (41) 

This relation is also motivated by the duality |T3]. The results of our cal- 
culations of inelastic cross section af^ e at ^fs = 1 TeV based on formulae 
(}40j) , (jUj) are presented in the second row of Table [TJ 

Table 1: Cross sections of the processes induced by graviton exchanges in i-channel 
(second row) and s-channels (third row) at ^/s = ITeV for different numbers of 
extra dimensions d (in pbarn). 



d 


2 


3 


4 


5 


6 


e + e~ -> e + e~ + X 


1.06 -10 3 


1.10 -10 2 


1.78 -10 1 


3.84 


1.02 


e+e" - // 


9.3 


3.7 


2.0 


1.3 


0.9 



These cross sections are larger than the cross sections of the processes 
induced by massive graviton exchanges in s-channel (at least, for d ^ 6). 2 For 
defmiteness, consider matrix element for fermion pair production e + e~ — > ff: 



m = g n t; u p^tI, Yl — T-jr- ( 42 ) 

nf+...n2>0 8 - J2 -A 

i=i R c 

Here P^ ua ^ is a tensor part of a graviton propagator, while Tfjj/' is the energy- 
momentum tensor of field e(f) f7JE3- The sum in (}4*2*)) diverges for d ^ 2. 
It can be estimated if one convert it into an integral and introduce an explicit 
ultraviolet cut-off M s . Then we get for d > 2: 



H 9 



n? + ...n^0 S-J2 T y 

i=i a. 



(rf-2)r(d/2)(4vr) d / 2 

(43) 

Rj 

r(l + rf/2)(4vr) d / 2 s ' V s 



In hadronic physics, the phcnomcnological parameter so ~ l/a'(0), where a'(0) ~ 1 
GeV -2 is the slope of hadronic Regge trajectories |12j . 

2 It is worth to note that, generally, QFTs for d > are not renormalizable. So, the 
following estimates are of illustrative character. 
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(an asymptotics at y^i <C M s was first found in [T%]). 

Taking into account that the sum in indices results in a factor proportional 
to s 2 , we obtain from (jUJ, 



~ As 2 — - x <^ 
\M D J 



1 

W 
1 



(44) 



where A = 0(1) has opposite sign for small and large y^s. Two asymptotics 
(JUJ) are well-matched at ~ M s . Thus, at y/s > M s we arrive at the 
expression 



2 N r fM a \W s 



where N c represents the number of colors of the final state. The result of 
numerical calculations by using formula (J43|) is presented in the third row of 
Tabled 

To compare, hadronic SM background in e + e~ annihilation (e + e~ — > 
e + e~+ hadrons), including the effects due to the (anti)tagging of the electron 
and accounting for all available data on 77 collisions, was estimated to be 

af a f(V~s = ITeV) ~ (2.7-4.0) • 10 4 pb. (46) 
The SM processes with different final states (Ylq^tt <?<?, W + W~, ti, x + X~ '> 

/i^/i^, Zh, etc.) have cross sections which are less than 1 pb at y/s = 1 
TeV (see, for instance, Fig. 1.3.1 from Ref. [16J). The highest rate has the 
process e + e~ — > J2 q ^t ( l ( l^ ^ s cross section is about 0.7 pb. 

Our goal is to find the process in which gravity forces can dominate SM 
interactions. Such a process has to obey the following requirements: (i) 
colliding energy is much larger than ~ 1 TeV, (ii) SM cross section 
does not rise rapidly in s. The best candidate is the scattering of ultra- 
high energy (UHE) neutrinos off the nucleons. These neutrinos is a part of 
ultra-high energy cosmic rays (UHECR) with energy E > 10 18 eV, which 
are dominated by extragalactic sources of protons [Hj. It is a detection 
of UHE neutrinos that can help us to discriminate between different origin 
of UHECR. For instance, in cosmological ( "bottom- up" ) scenarios, neutrino 
fluxes are almost equal to gamma ray fluxes. In astrophysical (acceleration) 
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approach, the neutrino flux is only a fraction of the gamma ray flux and is 
modified due to a propagation of cosmic rays before they reach the Earth. 

The cosmic neutrinos with extremely high energies E > 10 20 eV are 
also believed to explain so-called Greisen-Zatsepin-Kuzmin (GZK) cut-off of 
UHECR spectrum [TSj (see below). During UHECR propagation, the protons 
scatter off cosmic microwave background (CMB): 

P + Icmb^N + tt. (47) 

Taking into account that a typical CMB photon energies are 10 -3 eV, one 
can obtain that the nucleon interaction length drops to about 6 Mpc at 
GZK bound of E GZK ~ 5 • 10 19 eV [H2|. 3 The observation of UHECR at 
E > Eqzk is a serious problem for theories in which the origin of CR is 
based on acceleration of charges particles in astrophysical objects. Due to 
the energy losses (say, through process (jUj)), the UHECR particles cannot 
originate at distances larger than 60 Mpc from the Earth. On the other 
hand, all potential astrophysical sources of UHECR events are far beyond 
this distance. 

At the same time the process ()47|) is an origin of so-called cosmogenic 
neutrinos due to a consequence decay of charged pion as 7r ± — > /i ± z/ M , 
/i 1 * 1 — > v e v^. The fraction of the proton energy carried by the neu- 
trino is E v /E p pa 0.05 and independent of E p . The cosmogenic neutrino 
flux was first estimated in [21 j . [2*2*] . More recent estimates can be found 
in Refs. [221-1221 ■ The predicted fluxes depend on the evolution parame- 
ter m and on the value of the redshift z, and lye in the range: El <$> v ~ 
(0.5 ■ 10~ 9 - 10~ 8 ) GeV cm" 2 s" 1 sr" 1 at E v = 10 20 eV {u = fy, z/ e ). 4 

There are, however, other possible origin of UHE neutrinos. It is usu- 
ally anticipated that ~ $p M ~ $ y . We present below the total flux 
of muonic neutrinos and antineutrinos in a number of models at E v = 10 20 
eV. In the active galactic nuclei (AGN), the dominant mechanism for neu- 
trino creation is the accelerated proton energy loss due to pp or p7 in- 
teractions |2*7j . Note, AGN produce a large fraction of the gamma rays 
in the Universe, and their spectra agree with the prediction that gamma 
rays are produced by hadrons. In the AGN approach it was obtained that 

3 Below Eqzk, the dominant energy loss for the proton is due to the process pj CMB — > 
pe + e~, down to the threshold energy of 4.8 • 10 17 eV. 

4 Some cosmic ray protons with energies above 10 20 eV are converted into neutrons by 
pion photo-production. The neutrons decay again into protons during their propagation 
producing electronic anti-neutrinos. This mechanism is important at E Ve < 10 17 eV. 
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El $„ ~ 0.3 • 10" 8 GeV cm -2 s" 1 sr" 1 j27j. In Z-busts scenario, cosmic neu- 
trinos with extremely high energies [E v > 4 • 10 21 (leV/m u ) eV) collide with 
relic neutrinos 28J. If the masses of the background neutrinos m u are of 
several eV, the cosmic neutrinos initiate high energy particle cascades which 
can contribute 10% to the observed cosmic ray flux at energies above the 
GKZ cut-off (one of the main processes is a resonant uu collision via Z- 
bozon). The neutrino flux is El ~ 0.3 ■ 10" 6 GeV cm" 2 s" 1 sr" 1 ["""""j- In 
the so-called topological defect models [21"" , UHECR are produced via decays 
of supermassive X-particles related to a grand unification theory. The ex- 
pected neutrino flux is about El § u ~ 0.5 • 10" 6 GeV cm" 2 s" 1 sr" 1 [21]. In 
gamma-ray busts (GRB) model jHUj, the neutrino flux is strongly suppressed 
at E u > 10 19 eV, since the protons are not expected to be accelerated to 
energies much larger than 10 20 eV. 

It is worth to mention model-independent upper bounds on the intensity 
of high energy neutrinos produced by photo-meson interactions. If the size 
of cosmic ray source is not larger than photo-meson free path, the upper 
limit is (for evolving sources) 4.5 ■ 10" 8 GeV cm -2 s" 1 sr" 1 [3T]. However, 
for optically thick pion photoproduction sources, the upper limit is less re- 
strictive: 2.5 • 10" 6 GeV cm" 2 s" 1 sr" 1 [321 . Note, the considerably higher 
flux of cosmogenic neutrinos was obtained in Ref. [33]. The cosmogenic flux 
is the most reliable, as it relies only on two assumptions: (i) the observed 
extremely high energy cosmic rays (EHECR) contain nucleons, (ii) EHECR 
are primarily extragalactic in origin. 

One possible way to resolve the GZK puzzle 5 is to assume that the pri- 
mary UHECR particles are neutrinos which deposit a part of their energy 
to proton fragments in i/iV-interactions. Unfortunately, the SM neutrino- 
nucleon cross sections are not large enough to resolve the problem. Indeed, 
at 10 16 eV < E v < 10 21 eV the conventional contributions from charged and 
neutral current z/iV-scattering can be parameterized by [3*4] 



5 Note, however, recent paper [201, in which it is argued that the data from Fly's Eye, 
HiRes and Yakutsk cosmic ray experiments are consistent with the expected suppression 
of cosmic ray spectrum above 5 • 10 19 eV. The AGASA data show an excess in this region. 




(48) 
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The total SM cross section for z/iV-scattering has practically the same mag- 
nitude and energy dependence at energies under consideration ["""""] . Putting, 
E v = 10 21 eV in (jig)) , we get an estimate crg^ ~ 3.55 • 10 5 pb. Such a value 
of neutrino-nucleon cross section is two small to be relevant to the GZK 
problem. 

So, interactions beyond the SM 6 are needed in order to explain possible 
excess of the UHECR flux. One possibility is high-energy scattering mediated 
by gravitational forces in theories with compact extra dimensions [3l"]-j39j. 
In a number of papers [!""""]- [IE] it was shown that in a model with extra 
dimensions the neutrino-nucleon cross section can be enhanced by a black 
hole production. The corresponding cross sections were estimated to be one 
order of magnitude or more above the SM predictions (gHJ) at E v > 10 18 eV. 

In papers |40j . [4*2] the opportunities were considered to search for black 
hole signatures by using neutrino telescopes such as AMANDA/IceCube, 
Baikal, ANTARES or NESTOR. The expected black hole production cross 
section is around 10 6 pb for M^ff = Md = 1 TeV, where Mj"*jp is a minimal 
mass of produced black hole. 7 

Another possibility, which we will concentrate on, is an observation of 
air showers triggered by UHE neutrino interactions. The technique used for 
studying extensive air showers of UHECRs or UHE neutrino is the detection 
of shower particles by ground detectors, or the detection of fluorescence light 
produced by the shower. The first technique is used by one of the largest 
operating AG AS A experiment, while the second one was developed for Fly's 
Eye (HiRes) detector. The largest project under construction is the Pierre 
Auger Observatory [|7] . It will consist of two sites, each having 1600 particle 
detectors overlooked by four fluorescence detectors. For a detail study of 
extensive air showers with energy above 10 18 eV, 10% of the events will be 
detected by both ground array and fluorescence detectors. 

It is worth also mentioned space-based experiments EUSO and OWL 
which will be sensitive to CRs with energies above 10 19 eV. The future of the 
neutrino astronomy may be related with radio frequency detectors, such as 
RICE and ANITA. 

The neutrino interaction length is given by (in units of km water equiva- 

6 Thcrc is, however, a possibility that SM instanton-induced processes may give a large 
neutrino-nucleon cross-section [*"""| • 

7 The production rate of black holes depends on the number of extra dimensions and, 
essentially, on the ratio M^^ 1 /Mo- 
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lent, 1 km we = 10 5 g cm 2 ) 



L V {E V ) ~ 1.7-10 
For typical black hole production cross section ct^n — 

10 6 pb, we get L v = 17 

km we. This interaction length is much larger than the vertical Earth's 
atmospheric depth, which is equal to 0.01 km we. The atmospheric depth 
for neutrinos transverse (almost) horizontally is 36 times larger. That is why, 
it was proposed to search for uniformly produced quasi-horizontal showers 
at ground level [4*5]. 8 

The Fly's Eye and AGASA Collaborations have searched for deeply pene- 
trating quasi-horizontal air showers, with the depth L s h > 2500 g cm -2 . The 
probability of cosmic protons and gamma rays initiating air showers deeper 
than 2500 g cm~ 2 is about 10 -9 . Thus, any shower starting that deep in the 
atmosphere is a nice candidate for a neutrino event. 

The non-observation of such events puts an upper limit on the product of 
the neutrino differential flux, $^ = (l/4ir)dN v /dE v , times neutrino-nucleon 
cross section. The Fly's Eye Collaboration gives the bound which can be 
parametrized by 

/ F \ -1.48 

t$ v a vN ){E v ) < 3.74- 10- 42 ( ^-j GeV" 1 s" 1 sr" 1 , (50) 

while the upper limit from Ref. |SZ| can be recast as follows 

(*„ a vN ){E v ) < 10- 41 y- 1 ' 2 15 GeV- 1 s" 1 sr" 1 , (51) 

where y is an average fraction of the neutrino's energy deposited into the 
shower. The inequalities are valid in the range 10 8 GeV < E v < 10 11 GeV, 
provided a u jq{E v ) < 10 fib. 

Let us now estimate the neutrino-nucleon cross section in our approach. 
The neutrino scatters off the quarks and gluons distributed inside the nucleon 

8 At large zenith angles a background from hadronic cosmic rays is negligible, since the 
showers initiated by hadrons are high in the atmosphere due to very short interaction 
length of the proton. Around 10 20 eV, the hadronic mean free path is only 40 g cm~ 2 , 
and gamma-rays of such energy have interactions lengths of 45-60 g cm -2 . 



lpb 



o~vn{E v 



km we. 



(49) 
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(see a comment after formula (|53|) h Then the cross section is presented by 

1 

o£(a)= I dxJ2fi&» 2 >in(s), (52) 

J i 

where /j(x,/i 2 ) is a distribution of parton i in momentum fraction x, and 
s = xs is an invariant energy of a partonic subprocess. In our approach, 
the partonic cross section a in (s) is defined via the eikonal ()4U|1. As it follows 
from (|4"U|). x(s, b) is small at s < M|,, and we can put x m i n = M^/s in (J35J). 
At s ^ Af£, the main contribution comes from the region: 

b 2 < b 2 max (VI) = 4*2,(5) [Hs/Ml) + 1]. (53) 

We choose the neutrino energy E u to be 10 17 eV, 10 18 eV, 10 19 eV, 10 20 eV, 
and 10 21 eV. The invariant energy of uN collision is then 14 TeV, 43 TeV, 
137 TeV, 433 TeV and 1370 TeV, respectively Since b max (^s = 1370 TeV) ~ 
3-10~ 2 GeV -1 = 6-10 -3 fm (for 2 ^ d ^ 6), our assumption that the neutrino 
interacts with the proton constituents is well justified. 

We use the set of parton distribution functions (PDFs) from paper [201 
based on an analysis of existing deep inelastic data in the next-to-leading 
order QCD approximation in the fixed-flavor-number scheme. The extraction 
of the PDFs is performed simultaneously with the value of the strong coupling 
and high-twist contributions to structure functions. The PDFs are available 
in the region 10~ 7 < x < 1, 2.5 GeV 2 < Q 2 < 5.6 • 10 7 GeV 2 jSOJ - We 
take the mass scale in PDFs to be /i — l/b max (\/§), with b max defined by 
equation The result of our calculations of a^{E v ) is presented in 

Table |21 9 These neutrino-nucleon cross sections do not violate experimental 
upper bounds (|50|). ([ST]). 

Note, the total SM cross sections for (z/ + ^-scattering defined by formula 
IP are equal to 6.27 -10 3 pb, 1.45 -10 4 pb, 3.35 -10 4 pb, 7.72 -10 4 pb, and 
1.78 -10 5 pb, respectively. Thus, the SM interactions become comparable 
with (larger than) gravity interaction for d = 3 -=-4 (for d ^ 4 H-5), depending 
on the neutrino energy E u . 

The number of horizontal hadronic air showers with the energy E s h larger 
than a threshold energy E t h, initiated by the neutrino-nucleon interactions, 

9 The SM contributions to neutrino-nucleon cross sections are not included in Tabled 
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Table 2: Inelastic neutrino- nucleon cross section for the graviton induced scatter- 
ing at fixed neutrino energy, E u , for different numbers of extra dimensions d (in 
pbarn). 



d 


2 


3 


4 


5 


6 


E v = 10 17 eV 


8.63 -10 4 


5.63 -10 3 


5.53 -10 2 


7.16 -10 1 


1.13 -10 1 


E v = 10 18 eV 


6.53 -10 5 


3.39 -10 4 


2.47 -10 3 


2.26 -10 2 


2.41 -10 1 


E v = 10 19 eV 


4.20 -10 6 


2.05 -10 5 


1.21 -10 4 


8.59 -10 2 


6.99 -10 1 


E v = 10 20 eV 


2.05 -10 7 


1.32 -10 6 


7.06 -10 4 


4.29 -10 3 


2.94 -10 2 


E v = 10 21 eV 


8.74 -10 7 


7.47 -lO 6 


4.56 -10 5 


2.52 -10 4 


1.52 -10 3 



is given by 

N s h(E sh ^ E th ) 



TN A 



dE v § V {E V ) a 9 ™ v (E v ) A(E V ) 6(E U - E, 



th j 



+ J2 j dE ^ $ ^(^) <%(E Ut ) A{ Vl E Vi ) Q{ Vl E Vi - E th ) 



(54) 



=e,/j.,T 



where Na = 6.022 • 10 23 g _1 , T is a time interval, and A is a detector accep- 
tance (in units of km 3 steradian water equivalent). The quantity ^ Vi {E Vj ) in 
J5U) is a flux of the neutrino of type i, and <& V {E V ) = J2i= e « T ^^( E ^)- W 
The inelasticity yi defines a fraction of the neutrino energy deposited into 
the shower in the corresponding SM process (see below). 

The AGASA acceptance for deeply penetrating quasi-horizontal air show- 
ers with zenith angles 9 > 75° can be found in the second paper of Refs. [33] . 
It rises linearly in E s h in the interval 10 7 GeV < E s h < 10 10 GeV, while 
in the ultrahigh high-energy region the acceptance is constant and equal to 
A(E sh ^ 10 10 GeV) w 1.0 km 3 we sr ggj. 

The neutrino acceptance of the Pierre Auger detector is roughly 30 times 
larger, taking into account the ratio between Auger and AGASA surface 
areas. The acceptance of the Auger ground surface array has been studied in 
details in Ref . jM] , while the acceptance of fluorescence detector to neutrino- 
like air showers with the large zenith angles was calculated in Refs. 
[5*4] . The Auger observatory efficiency is high, since the low target density 
in the atmosphere is compensated by the very large surface area of the array 



^Both neutrino and antineutrino are everywhere included in the sum. 



17 



(each side of it covers an area of 3000 km 2 ). The highest efficiency for quasi- 
horizontal shower detection is expected at E sh > 10 9 GeV |q"T] . 

The number of extensive quasi-horizontal showers induced by the neutri- 
nos with energy larger than some threshold energy Eth, which can be detected 
by the array of the southern site of the Pierre Auger observatory, is presented 
in Table El The neutrino-nucleon cross section a 9 J£ v in Eq. (J53j) describes 
the contributions from the reggeized KK gravitons. The cosmogenic neutrino 
flux is from Refs. assuming a maximum energy of E max = 10 21 eV for 
the UHECR. The acceptance of the Auger detector is taken from Ref. [ST] 
(it is not assumed that shower axis falls certainly in the array). 

Table 3: Yearly event rates for nearly horizontal neutrino induced showers with 
^zenith > 70° and E s h E t h for cosmogenic neutrino flux from [2Hj at three values 
of threshold energy Eth- Number of events corresponds to one side of the Auger 
ground array. 



d 


2 


3 


4 


5 


6 


E th = 10 8 GeV 


34.88 


2.00 


0.32 


0.21 


0.20 


E th = 10 9 GeV 


30.21 


1.66 


0.21 


0.12 


0.12 


E th = 10 10 GeV 


13.16 


0.74 


0.062 


0.025 


0.022 



The neutrino-nucleon inelastic interactions induced by gravireggeons re- 
mind the SM neutral currents events. We assume that such events should 
result in hadronic dominated showers without leading lepton. That is why, 
we choose the inelasticity to be equal to unity 11 for the events induced by 
gravireggeon exchange (the first term in the RHS of Eq. (|54|) ) . We have 
also put y e = 1 for the SM charged current interactions initiated by elec- 
tronic neutrino, while for the SM neutral interactions initiated by u e and for 
v^/v r -e vents we have taken y e = y^ = y T = 0.24, following the calculations 
presented in Ref. [52] . 

The neutrino event rates are expected to be much higher for the neutrino 
fluxes obtained in "optimistic" scenarios considered in Refs. As an 

example, we have presented the yearly event rates for the Z-burst scenario in 

11 The estimates from Ref. are not applicable in our case, since in [^j an energy 
loss in elastic neutrino-nucleon cross section induced by "bare" gravitons was considered, 
while we deal with inelastic cross section in the gravireggeon model. 



18 



Table. 4. One can see from Tabled that the main contribution to the shower 
rate comes from the region of extremely high neutrino energies [E v > 10 10 
GeV). It can be understood as follows: at UHEs, the neutrino flux times E v 
varies slowly in E v in the Z-burst model (up to 2.5 • 10 12 GeV), while both 
the acceptance of the Auger array and "gravitational" part of the neutrino- 
nucleon cross section rise with the neutrino energy (see Table |2J) . 

Table 4: The same as in Table |3 but for the Z-burst neutrino flux from |28j . 



d 


2 


3 


4 


5 


6 


E th = W GeV 


12.60 ■ 10 2 


11.53 ■ 10 1 


9.26 


1.90 


1.50 


E th = W GeV 


12.59 • 10 2 


11.53 • 10 1 


9.26 


1.90 


1.50 


E t h = 10 10 GeV 


12.55 • 10 2 


11.51 • 10 1 


9.20 


1.85 


1.44 



The calculations of the yearly event rates in the energy interval 10 8 GeV ^ 
E sh ^ 10 11 GeV in the Z-burst scenario result in 44, 2.7, 0.38, 0.26, and 0.25 
for d = 2, 3, 4, 5 and 6, respectively. Remembering that combine results from 
AGASA and Fly's Eye imply an upper bound of 3.5 at 90% CL from quasi- 
horizontal neutrino events [12] , an d taking into account that the AGASA 
acceptance is roughly 30 times smaller than the Auger acceptance, we can 
conclude that the Z-burst neutrinos do not violate bounds (}50|) . (|5*T|) in our 
scheme for d ^ 3. 12 



3 Conclusions 

In the model with compact extra spacial dimensions, we have calculated 
the contribution of the KK gravireggeons into the inelastic cross section of 
the high energy scattering of both D-dimensional and four-dimensional SM 
particles. The usually adopted summing non-reggeized gravitons leads to 
a divergent sum in KK-number n (for D ^ 6). In our approach, on the 
contrary, the contribution of gravireggeon with the KK-number n to the 
eikonal is exponentially suppressed at large n. As a result, the corresponding 
sum in n is finite, and it can be analytically calculated. 

12 We do not discuss here cosmological bounds on the number of extra dimensions pQ 
(see also jSS] and references therein). 
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In the case when the SM fields propagate in all D dimensions, the de- 
pendence of the inelastic cross section on invariant energy y/s appeared to 
be similar to the upper limit for the total cross section obtained previously 
for the SM in the D-dimensional flat space-time without gravity. When, on 
the contrary, only gravity lives in extra dimensions, the imaginary part of 
the eikonal is derived in a closed form, which depends (except for y/s and 
the impact parameter b) on the number of extra dimensions d = D — 4 and 
their size R c . 

We have estimated the event rate for the quasi-horizontal air showers, 
induced by the interactions of UHE neutrinos with nucleons, which can be 
yearly detected by the ground array of the Pierre Auger observatory. It 
decreases rapidly if d varies from 2 to 5. For d — 4, we expect 10 events per 
year for the neutrino flux predicted in the Z-burst model. For the cosmogenic 
neutrino flux, gravireggeon induced interactions do not increase the event rate 
significantly with respect to the number of the neutrino events calculated in 
the SM, except for the case d < 3, which is likely to be excluded by the 
cosmological data. 
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